Intrinsic Origin of Negative Fixed Charge in Wet Oxidation for Silicon Carbide 
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We demonstrate on the basis of first-principles calculations that the formation of carbonate-like 
moiety in Si02 could be the intrinsic origin of negative fixed charge in SiC thermal oxidation. We 
find that two possible origins for the negative fixed charges are 0-lone-pair state and a negatively 
charged CO3 ion in Si02. Such CO3 ion is able to be formed as a result of the existence of residual 
C atoms in Si02, which are expected to be emitted from the interface between SiC and Si02, and 
the incorporation of H atoms during wet oxidation. 



SiHcon carbide (SiC) has recently attracted a great deal 
of attention as a promising material for the next genera- 
tion of high-power and high-frequency electronic devices. 
One of advantages of SiC for practical application is the 
formation of insulating silicon oxide (Si02) layers on SiC 
by thermal oxidation, being important for the fabrica- 
tion of SiC-based metal-oxide-semiconductor field-effect 
transistors (MOSFETs) [1-3]. There are two types of 
oxidation processes; dry and wet oxidation. The dry oxi- 
dation process of SiC is performed in O2 atmosphere. On 
the other hand, the wet oxidation process is conducted in 
O2 and H2O atmosphere, and, as compared to the case 
of dry oxidation, the wet oxidation process is known to 
effectively reduce the deep donor-type interface states for 
the p-type MOS capacitors, leading to the improvement 
of the channel mobility of the MOSFET 0, |5|. Further- 
more, wet oxidation yields high oxidation speed, twenty 
five times larger than that of dry oxidation, being also 
advantageous for realistic applications toward mass pro- 
ductions [6]. 

It is, however, serious issue encountered in wet oxida- 
tion. This treatment leads to surprisingly large amount 
of negative fixed charges near the SiC/Si02 interfaces 
[5.,^]. The wet-oxidation-induced interface traps between 
SiC and thermal oxide Si02 have been investigated exper- 
imentally by high-frequency C-V measurement of MOS 
capacitors [H, [H . It has been reported that fixed charges 
near Si02/SiC interfaces give large influences on thresh- 
old voltage and cannel mobility of SiC-MOSFET. No- 
tably, it has been pointed out that in wet oxidation the 
interface fixed charges increase up to —13x10^^ cm~^ 
(wet reoxidation annealed after wet oxidation) in p-4H- 
SiC MOS capacitor [5]. Non- negligible negative fixed 
charge also appears during dry oxidation [9', ^10^ . Nev- 
ertheless, despite the extensive experimental outcomes, 
the origin of the negative fixed charge remains unclear 
for the last decade. 

In this letter, we propose to clarify the atomistic ori- 
gin of negative fixed charge appearing in SiC thermal 
oxidation processes on the basis of first-principles cal- 
culations in the framework of density functional theory 
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(DFT). We here focus on residual C atoms in Si02 which 
are expected to be emitted from the interface between the 
SiC and Si02. The existence of such residual C atoms 
is supported through a comparison with Si thermal oxi- 
dation. Indeed, it has been reported experimentally and 
theoretically that one-third Silicon atoms are inevitably 
emitted from the interface to release the stress induced 
during Silicon oxide growth [11I5 [HI • 

Our calculations demonstrate that residual C atoms in 
Si02 leads to the formation of carbonate-like ion assisted 
by incorporated H atoms. We find that the existence of 
the carbonate-like ion stabilizes 1— charge state, being 
the intrinsic origin of negative fixed charge in Si02/SiC 
interface. 

We considered three types of structural models based 
on 72-atoms a-quarts Si02. The first model was con- 
structed by substituting one C atom for Si, while the 
second and third models involve one and three H atoms, 
respectively. The lattice vectors were set to be ai = 
(1, — >/3, 0)a, a2 = (1, >/3, 0)a, and as = (0, 0, 2)c, where 
a = 4.944 A and c = 5.363 A, which was optimized for 
Si02 bulk in the a-quarts phase. All calculations were 
performed using DFT with the Perdew-Burke-Ernzerhof 
generalized gradient approximation p!3l-[l5| . We used ul- 
trasoft pseudopotentials for Si, O and C and normcon- 
serving pseudopotentials for H [161, [13 • We used a plane- 
wave basis set with a cutoff energy of 36 Ry. We sam- 
pled 2x2x2 /c-points in the Brillouin zone (BZ) for the 
BZ integration; using 3x3x3 /c-points in the BZ gives a 
difference of 0.01 eV in total-energy for the system, assur- 
ing required accuracy in the present calculations. After 
the optimizations, all the atomic forces were less than 
0.05 eV/A. We used the background ion charge correc- 
tion proposed by Blochl to calculate the charged state 

Figure [1] shows the the formation energy diagram of 
the models including one C atom, one C atom and one H 
atom, and one C atom and three H atoms, respectively. 
We found that a defect with only C atom takes the neu- 
tral charge state for the whole range of the Fermi energy 
(the upper panel in Fig. [U (a)). However, the system 
including H atoms takes the 1— charge state for nearly 
the whole range of the Fermi energy (the upper panels 
in Fig. [U (b) and (c)). Particularly, if we assume that 
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FIG. 1. (Color online) Formation energy diagram and optimized structures for the models including (a) one C atom, (b) one C 
atom and one H atom, and (c) one C atom and three H atoms in the unit cell. Dashed, solid, and dotted lines in the formation 
energy diagram show neutral, 1 — , and 2— charge state, respectively. Optimized structures are shown for neutral state in the 
panel (a), and for 1— charge state for the panels (b) and (c). Distances are measured in units of A. 



the valence band offset parameter between 4H-SiC and 
Si02 is 2.9 eV [19], the 1— charge state appears as the 
most stable state for all over the range of SiC band gap 
region. These results indicate that negative fixed charge 
is formed when both of C and H atoms exist. 

Notably, we found that the structural origin of the 
stabilization for 1— charge state is the formation of 
carbonate-like ion. The lower panels in Fig. [1] shows 
the optimized structures for each model. It is clearly 
shown in the figure that carbonate-like (COs-like) moi- 
ety is formed in Si02 in the H-atom-included cases (Fig. 
[1] (b) and (c)). The formation of COs-like configuration 
is ascribed to H atoms that can cut Si— O bonds near the 
substitutional C atom, leading to the formation of strong 
sp^-like C-0 bonds. These results provide clear indica- 
tion that C atoms originating from Si02/SiC interface 
leads to the formation of sp^ network in Si02, which is 
assisted by H atoms. 

There are two main factors for the stabilization of the 
COs-like moiety. The first one is partial resonance in the 
COs-like configuration. Indeed, as shown in the lower 
panels of Fig. [1] (b) and (c), its C-0 bonds are charac- 
terized as partially double bond and fully double bond, 
as compared with the typical lengths of the C-0 single 
and double bonds (1.42 and 1.21 A, respectively). The 
second factor is the formation of hydrogen-bonding be- 
tween the COs-like part and the nearby Si-O-H moiety 
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FIG. 2. Schematic illustration of structural change for COs- 
like moiety. 



(Fig. 111(c)). 

The detailed structure of the COs-like moiety depends 
on the number of incorporated H atoms in the system. 
In the case of the one H incorporation, two of three O 
atoms in the COs-like part form covalent bonds with Si 
atoms (the lower panel in Fig. [1] (b)). The COs-like 
structure has C-0 bonds of 1.26, 1.36, and 1.36 A, and 
0-C-O angles of 113, 122, and 125 degrees. On the other 
hand, in the case of the incorporation of three H atoms, 
only one O atom in the COs-like moiety form a covalent 
bond with Si atom (the lower panel in Fig. [1] (c)). In 
this case, C-0 bonds become 1.27, 1.29, and 1.36 A, and 
0-C-O angles of 116, 121, and 123 degrees. These struc- 
tural change of COs-like moiety can be understood by 
considerations about its local bond network (Fig. [2j). As 
compared with an isolated COs-like ion, Si-0 bond for- 





FIG. 3. (Color online) Energy levels around energy gap at 
the r point and the distribution of wave functions for the 
system including one C atom and one H atom in 1— charge 
state. Energies are measured from the Si02 valence band 
top. The isovalues for the wave- function distribution are + 
0.12 (gray) and - 0.12 (black) (e/A^)^/^ 



mation leads to a structural change of COa-like moiety, 
depending on the number of Si-0 bonds. This chemical 
picture also strongly suggests that COa-like part is in 1 — 
charge state when one Si-0 bond is formed, while it is in 
neutral charge state when two Si-0 bonds are created. 

Indeed, a detail analysis of the electronic structure re- 
veals that there are two origins for negative fixed charge, 
depending on the number of incorporated H atoms in the 
system. Figure [3] shows the energy level and the distri- 
bution of wave functions for the system including only 
one H atom in 1— charge state. There are five occupied 
gap states named a, /3, 7, and e in the Si02 band gap 
region. The energy levels for the a, P and 7 states ap- 
pear in the range of ~ 1 eV from the Si02 valence band 
top. These states are distributed primarily on COa-like 
moiety and the nearby O atom. On the other hands, 
the S and e states emerge above ~ 2 eV from the Si02 
valence band top. These states are distributed only on 
O atom with lone-pair characteristics. These results in- 
dicate clearly that COa-like moiety is in neutral charge 
state, and an extra electron is trapped into 0-atom lone- 
pair electron state, which leads to the intrinsic origin of 
negative fixed charge. 

In sharp contrast, the three- H- atom incorporation case 
exhibits completely different electronic structure. There 
are only two gap states appearing above ~ 1 eV from the 
Si02 valence band top. The corresponding wave func- 
tions are distributed strongly in COa-like moiety. This 
indicates clearly that an extra electron is trapped into 
the COs-like part, and the moiety is in 1— charge state. 
The COs-like ion is thus the intrinsic origin of negative 
fixed charge when enough H atoms are incorporated in 
the system, and this could be possible in wet oxidation. 
In particular, the energy level for this defect state is deep 
as compared with Si02 large gap. This stabilization is 
ascribed to the formation of hydrogen-bonding between 
the COs-like ion and the nearby Si-O-H moiety as men- 



FIG. 4. (Color online) Energy levels around energy gap 
at the r point and the distribution of wave functions for the 
system including one C atom and three H atoms in 1 — charge 
state. Energies are measured from the Si02 valence band top. 
The isovalues for the wave-function distribution are + 0.12 
(gray) and - 0.12 (black) (e/A^)^/^ 




FIG. 5. (Color online) Optimized structure for the system 
including one C atom and three H atoms in 1— charge state. 
Distances are measured in units of A. 



tioned above. These results strongly suggest that the de- 
fect state originated from COs-like ion does not become 
a charge trapping level, but plays a strong negative-fixed- 
charge appearing in wet oxidation. 

On the other hand, it has been experimentally reported 
that a few but non-negligible negative fixed charges ap- 
pear during dry oxidation [9^, lIO,]. In this case, the neg- 
ative fixed charge could be 0-lone-pair state. In reality, 
Si02 structure is amorphous. It is thus expected that a 
few O atoms which is not terminated by H can be formed 
near incorporated C atoms. This situation is similar to 
the case shown in Fig. [1] (b), where lone-pair electron 
state originated from non-terminated 0-atom appears in 
Si02 band gap (Fig. ^. 

Finally, we note charge-state-dependent behavior of H 
atoms near COs-like ion. As shown in Fig. [H (c), a 
H atom terminates Si dangling bond when COs-like ion 
is formed. This H atom is, however, transferred from 
Si atom to one O atom in the COs-like ion when the 
system is in neutral charge state (Fig. [5]), leading to the 
formation of Si dangling bond. 
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In conclusion, we have demonstrated that the forma- 
tion of carbonate-hke moiety in Si02 dielectrics could be 
the origin of negative fixed charges in the SiC thermal oxi- 
dation using first-principles calculations in the framework 
of DFT. Our calculations showed two possible atomistic 
origins for the negative fixed charges: 0-lone-pair state 
and a negatively charged CO3 ion. The CO3 ion is able to 
be formed as a result of the existence of residual C atoms 
in Si02, which are expected to be emitted from the in- 



terface between SiC and Si02, and the incorporation of 
H atoms during wet oxidation. 
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